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Abstract: We present an analysis and classification of vector and scalar fluctuations in 
a D3/D7 brane setup at finite termperature and baryon density. The system is dual to 
an = 2 supersymmetric Yang-Mills theory with SU{Nc) gauge group and Nj hyper- 
multiplets in the fundamental representation in the quenched approximation. We improve 
significantly over previous results on the quasinormal mode spectrum of D7 branes and 
stress their novel physical interpretation. Amongst our findings is a new purely imaginary 
scalar mode that becomes tachyonic at sufficiently low temperature and baryon density. 
We establish the existence of a critical density above which the scalar mode stays in the 
stable regime for all temperatures. In the vector sector we study the crossover from the 
hydrodynamic to the quasiparticle regime and find that it moves to shorter wavelengths for 
lower temperatures. At zero baryon density the quasinormal modes move toward distinct 
discrete attractor frequencies that depend on the momentum as we increase the temper- 
ature. At finite baryon density, however, the trajectories show a turning behavior such 
that for low temperature the quasinormal mode spectrum approaches the spectrum of the 
supersymmetric zero temperature normal modes. We interpret this as resolution of the 
singular quasinormal mode spectrum that appears at the limiting D7 brane embedding at 
vanishing baryon density. 
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1. Introduction and summary 

One of the successes of the AdS/CFT correspondence [1,2] and its generalizations is its 
application to the plasma phase of non-Abelian gauge theories [3] . This is of high interest 
because of its potential relevance for the description of the strongly-coupled quark-gluon 
plasma as it is created in Heavy Ion Collision at RHIC or in the near future at the LHC. 

The holographic modelling of the plasma phase invariably involves an asymptotically 
AdS black hole [4] . Of particular interest are the quasinormal modes of such black holes 
as they are mapped to the poles in the correlation function of the dual finite temperature 
field theory [5,6]. One aspect of this is a relation between the quasinormal frequencies and 
the hydrodynamic transport coefficients, e. g. the shear viscosity. 

An important generalization of the AdS / GET correspondence is the addition of flavor 
degrees of freedom in the fundamental representation of the gauge group. One convenient 
way of achieving this is via the addition of D7-brane probes to the ten-dimensional super- 
gravity background [7]. The meson spectrum can then be studied via the fluctuations of 
the probe brane [8,9]. 

Eor the configuration of a D7-brane probe added to the AdS Schwarzschild black hole 
background, a first order phase transition occurs between D7-brane probes either staying 
outside of the horizon or reaching down to it [10-12]. These two types of embeddings are 
called Minkowski or black hole embeddings, respectively. As discussed in [13], the first 
case corresponds to stable mesons, whereas in the second case the mesons are unstable. 
Eor stable mesons, the spectrum has been determined analytically at zero temperature 
in [8] . In the second case, the meson excitations may be identified with quasinormal modes 
and their finite width is related to an infalling energy boundary condition at the black 
hole horizon. The spectral functions for this configuration were first studied in [14]. If the 
temperature is high compared to the quark mass the quasinormal frequencies lie deep inside 
the lower complex half plane and the spectral function is smooth and without noticeable 
peaks. As the temperature is lowered the quasinormal modes move towards the real axis, 
producing rather well defined quasiparticle peaks in the spectral function. This behavior 
is particularly strong at finite baryon density. 

Eurthermore a tachyonic mode is present in the scalar sector close to the first or- 
der phase transition discussed above. Eor null momentum, the spectral function of this 
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configuration has been studied in [15] and the presence of the tachyonic mode has been 
argued for by way of a WKB analysis. We find indeed a so far overlooked purely imaginary 
mode in the scalar spectrum that crosses the real axis and therefore becomes unstable at 
a temperature that is in excellent numerical agreement with [15]. 

At finite baryon chemical potential, as obtained by considering a non-zero profile for 
the time component of the gauge field on the D7-brane, there is a rich phase structure for 
the D7-brane embedded in the AdS Schwarzschild background [16-20]. In these works the 
strict supergravity limit is considered, i.e. ol — > 0, and the results are thus perturbative on 
the gravity side. Non-perturbative effects such as worldsheet instantons have been studied 
in [21]. Taking these effects into account the first-order phase transition mentioned above 
becomes third order at a critical chemical potential [21]. In this paper however we will 
consider the strict supergravity limit only. 

The spectral functions for vector mesons at finite baryon density have been studied 
in [22] , where a particular turning behavior of the quasi-normal frequencies in the complex 
frequency plane was observed as function of the ratio of quark mass over temperature: At 
low values of Mg/T, the position of the quasinormal modes moves to smaller real parts 
when Mq/T increases, while for large values, it moves to larger real parts. For large Mq/T, 
the spectrum approaches the form of the supersymmetric case. The momentum dependence 
of the spectral function at finite baryon density has been investigated in [23,24]. 

Summary of results The purpose of the present paper is to present an in depth study 
of both scalar and vector modes at either finite baryon density or finite momentum. 

At zero density and zero momentum we determine the first and second quasinormal 
frequency of the scalar and vector fluctuations. We calculate their trajectories in the 
complex plane as we increase the quark mass over temperature ratio. In general these 
modes move along a curve where initially the real and imaginary parts can grow a little 
but then move continuously towards lower values. As shown in [25] the endpoints of 
these curves asymptote to a single point on the real frequency axis that is reached as the 
embedding reaches the limiting embedding. This happens when the brane just touches the 
black hole horizon^. There it has been argued that the spectrum of quasinormal modes 
becomes singular as they all coincide at one (real) energy value. It is however unclear if this 
limiting embedding can really be reached since the embeddings become locally unstable 
at rather high temperatures. Indeed in the scalar sector we find a mode that becomes 
tachyonic shortly after the system has become ther mo dynamically metastable due to the 
presence of the first order phase transition. Once the tachyonic mode is in the spectrum 
the embedding is not even metastable and so far it is not known what the true ground 
state would be or if such a state exists at all. 

At finite momentum but zero density, we find that at a critical momentum the system 
undergoes a crossover transition from the hydrodynamic behavior at long wavelengths to 

^Note that this hmiting embedding is often called critical embedding in the literature e. g. [12]. However 
in this paper we use the term 'critical embedding' for the embedding at which the phase transition occurs 
(cf. figure 1 a). 
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a collisionless behavior at small wavelengths. We observe this crossover in the longitudinal 
vector channel which has a hydrodynamic mode, i. e. a mode whose dispersion relation 
does not show a gap at zero momentum. This mode describes baryon charge diffusion and 
is thus called diffusion mode. The transition from the hydrodynamic to the collisionless 
regime is defined as the point where the imaginary part of the diffusion mode becomes 
larger than the imaginary part of the first quasinormal mode which is not a hydrodynamic 
mode. Additionally we observe that the first quasinormal mode of the scalar and vector 
fluctuations with increasing momentum show an increasing number of spirals in their tra- 
jectories parametrized by the quark mass over temperature ratio. The number n of these 
spirals apparently is correlated with certain "attractor" frequencies, i.e. the real part of 
the quasinormal modes asymptotes to a;„ G R as the D7-brane embedding approaches 
the critical embedding between the black hole and Minkowski embeddings. Although the 
asymptotic values a;„ lie deep in the unstable sector, the spirals mostly lie in the physical 
sector of the theory. So they are physical features of the (meta) stable theory. 

At finite quark density but zero momentum, we find the following: In the vector mode 
sector, there are two distinct movements of the quasinormal frequencies in dependence 
of the quark mass over temperature ratio. At small densities there is a turning point as 
already observed in [22] and as described above. This turning point is no longer present 
for large densities. The critical value for the normalized density is dc = 0.04. The relation 
between the normalized density d and the baryon density is given by 



where T is the temperature, Nf the number of flavors and A the 't Hooft coupling. In the 
scalar channel, the quasinormal frequencies show also distinct movements in dependence 
of the quark mass over temperature ratio. In this case we find three distinct behaviors 
separated by critical densities. As a function of Mq/T, the quasinormal modes perform 
either a left turn, a right turn or no turn in the complex frequency plane. We will study 
these different behaviors in detail. 

For the scalar channel, we find in addition that the tachyonic mode mentioned above 
is only present below a critical normalized density, d < 0.00315. We identify the unstable 
region in which the tachyon is present with an instability related to the first order phase 
transition found in the thermodynamics of the system. This behavior may be understood 
by comparing it to the first order phase transition of the well-known van-der-Waals gas, as 
described for instance in [26]. Figure 1 (b) shows the first order phase transition of the van- 
der-Waals gas, which follows the red line in the p-F-diagram as obtained from the Maxwell 
construction. Additionally there are two metastable branches: One of them reaches from 
the point where the phase transition begins when increasing the volume to the minimum 
denoted by the red plus sign in the figure. The second reaches from the maximum denoted 
by the red cross to the point where the phase transition ends when increasing the volume. 
The region between the two extrema is unstable since the pressure raises when the volume 
is increased. 



d 



2'2nB 





-4- 



Our calculations show a very similar behavior of the system considered here. We have 
calculated both the free energy as function of m oc Mq/T close to the phase transition, 
^ displayed schematically in Figure 1 (a), and the quasinormal modes for the same range 
of Mq/T. Both calculations show a similar structure of stable, metastable and unstable 
branches, with the numerical values for the boundaries of these branches agreeing up to 
four digits in both calculations. In particular in analogy to the van-der-Waals gas we find 
that in addition to the stable branches^ there are metastable branches close to the phase 
transition, denoted as overheated and undercooled^ in figure 1 (a). Furthermore we find 
an unstable branch which connects these two metastable branches. The metastable phases 
are stable against fluctuations, while on the unstable branch a tachyonic mode appears in 
the quasinormal spectrum. In particular we emphasize that at zero density the limiting 
embedding which has a conical singularity at the black hole horizon, and is often discussed 
in the literature, lies clearly deep inside the unstable region. Thus any observation found 
by using this embedding should be considered with great care since it does not correspond 
to a physical state. Furthermore we should stress that on the field theory side there exists 
a stable ground state for any combination of the mass m and chemical potential ji. 



Figure 1: (a) Sketch of the free energy F of the flavor fields versus the quark mass over temperature 
ratio m oc Mq/T close to the first order phase transition, (b) Pressure versus volume of the van-der- 
Waals gas. The red line marks the phase transition which is obtained by the Maxwell construction. 

We also present a qualitative analysis of the quasinormal spectrum at either finite 
momentum or finite density which uses the fact that the equation of motion for the fluctu- 
ations can be transformed into a Schrodinger equation. The Schrodinger potential analysis 
was introduced for instance in [13, 14,23,25]. 

A similar calculation of the free energy was performed in [16] . 

^Note that in [16] it was observed that for a small region of the stable phase in which the mass parameter 
m is just slightly larger than the critical one there is an charge instability. This instability is distinct from 
the instability discussed above. We do not observe this charge instability in our quasinormal mode analysis. 
This may be due to inadequacy of the our numerical methods at large mass over temperature ratios. 

*In the context of the D3/D7 system, the terms 'overheated' and 'undercooled' were introduced in [25]. 
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The results of this analysis are in qualitative agreement with the results obtained by 
stating the quasinormal modes listed above: In particular, at finite density we find a general 
feature in the Schrodinger potential: At low quark mass over temperature ratio there are 
only unbound scattering states which correspond to quasinormal modes. By increasing 
the quark mass over temperature ratio, a barrier forms in the Schrodinger potential [23] 
which leads to a local minimum. In [25] it was found for Minkowski embeddings that the 
Schrodinger potential is a box whose extent in the AdS radial direction coincides with the 
region filled by the D7-brane probe. Here we observe that for black hole embeddings at 
finite density, a barrier forms in the Schrodinger potential at the same radial position as the 
IR boundary of the box in the Minkowski case. This barrier separates the bulk of the AdS 
space from the horizon of the black hole. As the barrier increases when the ratio of quark 
mass over temperature m oc Mq/T raises, less energy can leak into the black hole. Thus 
'bound' states which correspond to 'normal' modes are formed. We study the appearance 
and the behavior of this barrier in detail. 

The paper is organized as follows: In section 2 we give a general introduction to 
the relation between quasinormal modes and hydrodynamics. In section 3 we introduce 
the D3/D7-brane setup. In the sections 4, 5 and 6 we determine the quasinormal modes 
of the scalar and vector fluctuations first at vanishing momentum and density and later 
at finite momentum or finite density. In every section we also calculate the corresponding 
Schrodinger potentials which we use to describe the qualitative behavior of the quasinormal 
modes. At the end of each section we summarize the physical results which we found. 

2. Quasinormal modes and Holography 

In this section we recall the definition of quasinormal modes of black holes and the role 
they play in determining the response of a holographic field theory close to equilibrium. 

Quasinormal modes of a black hole are distinct perturbations of the black hole solution. 
Roughly they can be understood as resonances of the black hole. However since the energy 
of the perturbation can leak into the black hole, these fluctuations are not normal modes 
and thus have been dubbed quasinormal. Their corresponding frequencies consist of a 
real and an imaginary part. As for the damped oscillator, the real part of the frequency 
essentially determines the energy of the fluctuations, while the imaginary part is responsible 
for the damping. In AdS spacetimes the quasinormal modes satisfy the following boundary 
conditions. At the horizon they are purely ingoing, whereas at the conformal AdS boundary 
they have an asymptotic behavior that corresponds to a normalizable mode. In this paper 
we determine the quasinormal spectrum of the D7-branes. The corresponding modes can 
be grouped in terms of their transformation properties under spatial rotations of SO{3). 
We consider scalar modes given by perturbations of the brane embedding as well as vector 
modes given by perturbations of the gauge field on the brane. As is well established by 
now, the quasinormal frequencies of the dual gravity theory can be identified with the poles 
of correlation functions in dual thermal gauge theories [5,6]. 

A system which is close to equilibrium can be described by linear response theory. 
There the effect of an external perturbation is given by the retarded two point function 



-6- 



folded against the source of the perturbation. By a Cauchy integration in the complex 
frequency plane the response can be written as a sum over the contributions of the different 
quasinormal modes (see e.g. [27]). Writing the time dependence as exp(— ■io't) we note that 
a relaxation towards equilibrium can only happen if all the quasinormal modes lie in the 
lower complex half plane. Following [27,28] the response caused by an external perturbation 
j{LO, k) can therefore be written as a sum over quasinormal frequencies 

k)) = i9{t) Rn JK, k) e-*^"*-^"* , (2.1) 

n 

where the quasinormal frequencies are ujn = On " i^n and their residues arc Rn- If a 
mode comes to lie in the upper half plane it results in an exponentially growing mode and 
therefore represents an instability of the system. As we will see in the following, such an 
instability does indeed occur in the scalar sector of the D7-brane fluctuations. 

Of particular interest is the hydrodynamic limit that considers only fluctuations with 
small frequency and large wavelength. This hydrodynamic expansion can be seen as an 
effective field theory where the degrees of freedom with large frequency and small wave- 
length are integrated out. In this effective field theory only the poles of correlation functions 
closest to the origin are important. More precisely, the modes with a dispersion relation 
obeying lim^_^Q u;(/c) = represent the hydrodynamic regime. Hydrodynamic transport 
coefficients such as the shear viscosity or the diffusion constants can be read off form these 
poles. Finally let us mention that recently it has been shown that the determinant of wave 
operators in some asymptotically black hole backgrounds can be written in terms of the 
quasinormal modes [29] . 

3. Holographic Setup 

We are interested in a large Nc thermal field theory, at finite baryon chemical potential 
including fundamental and adjoint matter. Our matter content is that of the A/" = 4 super 
Yang-Mills theory and a number Nj of A/" = 2 super Yang-Mills fundamental hypermulti- 
plets. We consider the quenched approximation with Nc ^ oo and Nf <^ Nc fixed. In this 
limit the theory stays conformal at leading order in Nc- 

The dual gravity setup is given by a stack of Nc D3-branes and Nf probe D7-branes. 
The metric generated by the D3-branes is a non-extremal AdS black hole background 
placing the field theory at finite temperature. The D7-branes allow for strings which have 
one end on a D3- and the other end on the D7-brane. These 3-7 strings correspond to 
quarks in the fundamental representation. The length of these strings, i.e. the separation 
between the stacks of D3- and D7-branes corresponds to the quark mass Mq in the field 
theory. On the world volume of the D7-branes we introduce a background gauge field A 
which generates a chemical potential ^ in the field theory at the boundary. Depending on 
the problem it is convenient to work in different coordinate systems. We now introduce 
the p-coordinates and the ^-coordinates. 
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3.1 The />-coordinate system 

We find it convenient to use the p-coordinates of [16] to write the AdS black hole background 
in Minkowski signature as 

ds' = 4-2 f-^di^ + f^^") + ( -V + , (3.1) 



2i?2 V / J \Q 

with dOg the metric of the unit 5-sphere and 

/(£>) = 1-4, /(f) = 1 + 4' (3-2) 

where R is the AdS radius, with 

R'^ = AirgsNc a'^ = 2A a'^ . (3.3) 

This type of radial coordinate is better suited for calculation of the quasinormal modes 
with the shooting method, where one places a cutoff at some large value of p. 

The temperature of the black hole given by (3.1) may be determined by demanding 
regularity of the Euclidean section. It is given by 

T=||. (3.4) 

In the following we may use the dimensionless coordinate p = g/ Qh, which covers the range 
from the event horizon at p = 1 to the boundary of the AdS space at p — > oo. 

To write down the DBI action for the D7-branes, we introduce spherical coordinates 
{r, $73} in the 4567-directions and polar coordinates {L,(j)} in the 89-directions [16]. The 
angle between these two spaces is denoted by (0 < < vr/2). The six-dimensional space 
in the 456789-directions is given by 



d£»^ + g'^dnl = dr^ + r'^d^l + dL^ + L'^dcf)'^ 

= dg^ + g^{d@^ + cos^ ed^^ ^ ^^^2 q^^^^ ^ 



(3-5) 



where r = ^sin0, g"^ = r'^ + and L = ^?cos0. Due to the symmetry, the embedding of 
the D7-branes only depends on the radial coordinate p. Defining x = cos 0, we parametrize 
the embedding by x = x{p) ^-iid choose (j) = ^ using the 0(2) symmetry in the 89-direction. 
The induced metric G on the D7-brane probes is then 

ds\G) = ^ (-£dt2 + /df^) + f ^^4^^^^^^' + - ^'^^^^ • (^-^^ 
The square root of the determinant of G is given by 



-G = ^g^ffil - x')\/l -X^ + g^idexy , (3.7) 



where is the determinant of the 3-sphere metric. 
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3.2 Brane setup and background fields in ^-coordinates 

Here we introduce the coordinates also used in [13]. This set of coordinates maps the 
compact interval [0, 1] to the distance between the conformal boundary and the black hole 
horizon. Such a compact radial coordinate is particularly well suited for the calculation 
of quasinormal modes using the relaxation method as explained in appendix B. We work 
in a non-extremal black-hole background generated by the stack of D3-branes giving the 
effective metric 

+ dni„ (3.8) 



_/(,.)dt2 + ^+df2 



m 

with f{z) = 1 — z^. Note that the horizon is located at z = 1, the AdS-boundary at 
z = 0. Frequencies and momenta measured in (3.8) are related to physical frequencies and 
momenta by {tOph, kph) = 7rT{uj, k), where T is the temperature. 
The D7-brane action reads 

S = j d^^y/-det{P[G] + 2TTa'F) , (3.9) 

where P[G] is the metric's pull-back on the Nf D7-branes. Here the induced metric is 
given by 



^ = -4? + - + ^dx- + sin^ e(z) dni 



(3.10) 



where Q{z) describes the D7-brane embedding in the AdS-Schwarzschild background and 
the embedding coordinates are = {t,x, z, 01,02,03). To second order in the field 
strength the D7-brane Lagrangian can be written in the following way 

£ = y/-P[G] \l + TT^a'^FabF"''] . (3.11) 

With the zero order term, we obtain the background equation of motion 

= 3cose{z)[-l + z^{-l + z^)e'{zf] (3.12) 
-zsin e(z)[(3 + z^)@'{z) + 2^2(1 _ /)(2 - z^)@'{zf - z{-l + z^)@"{z)] . 

The brane embedding can be found by integrating this equation from the horizon out to 
the boundary. As initial conditions one chooses xo = cos(0(l)) and demands regularity 
on the horizon. The quark mass can be read off from the behaviour at radial AdS infinity 
located at z = [13] 

Mg = ^my/XT, m = x{z = 0) . (3.13) 

We have plotted the mass parameter m as a function of xo = cos0(f3/f) in figure 2. 
The mass is not a single valued function of xo- Although m is the physical parameter of the 
brane embedding we prefer to quote xo values instead of m since much of our investigations 
will take place in the regime where m ceases to be single valued. The change in the sign of 
dm/dxo is also suggestive of an instability. Indeed the maximum of m = 1.31 is reached 
at Xo = 0.962, which is precisely the value from where an unstable mode appears in the 
scalar sector of the quasinormal mode spectrum. 
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Figure 2: Plot of the dimensionless mass parameter m vs the cosine xo = cos Go of the embedding 
angle 0(z) at the horizon, i.e. Oq = 0(z = 1). The mass is not a single valued function of xo on 
[0, 1]! It takes a maximum value of m = 1.31 at xo = 0.962. The horizontal line indicates the first 
order phase transition at xo = 0.939. 
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Figure 3: At finite charge density: The dependence of the dimensionless quark mass m = 
2Mq/\/\T on the hor izon value xo = linip^i x of the embedding as shown in [30]. 

Looking at figure 3 it is obvious that at finite charge density the unstable region 
disappears at a critical density near dc = 0.00315 and for larger densities the quark mass 
is a monotonously increasing function of the embedding parameter xo^ i-e- Xo ~ Mq/T for 
dc > 0.00315. This critical value will be confirmed in our analysis of the quasinormal mode 
spectrum below. 

The second order term produces the Maxwell equation 

d,{y^^P\G]F''') = 0, (3.14) 

where Fab = da Ah — dbAa- We can choose the Az = gauge and expand Aa in plane wave 
modes. Moreover if we write the equations in a gauge invariant way using the electric fields 
in longitudinal El = kphAo + oOp^Ai and transverse direction Et = LOphAj- the equations 
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of motion are 



E'L{z) + 



C{z) + 



/(z)(a;2-/(z)A;2) 



C{z) + 



E'l{z) + B{z){lj^ - fiz)k^)ELiz) = , (3.15) 
Et{z) + B{z){uj'^ - f{z)k^)ET{z) =0, (3.16) 



with 



C{z) = -^+2z{-2 + z*)&{zf 



Our results are summarized in the following three sections. 
4. Vanishing momentum and density 

Our analysis produced a considerable amount of data and we will not show all of it because 
of some redundancy in the results. Technical details on our numerical methods are also 
deferred to the appendices A and B. The idea for this and the following two sections 
is, to have example figures for each case of interest and a listing of all the effects we 
observe. We also provide a qualitative analysis by studying the correspondent Schrodinger 
equations. The discussion subsection in each of the cases is then devoted to the physically 
most interesting effects, i.e. the tachyon, diffusion mode, turning point. 

4.1 Transverse vectors 

The transverse equation of motion can be written in this simplified form 

E'^{z) + A^{z)E't{z) + i?(z)(u;2 - f{z)k^)ET{z) = , (4.1) 

where Ai{z) = C{z) + Close to the boundary of AdS (z — > 0"*"), the differential 

equation reduces to 

Et{z) - -^Et{z) = , (4.2) 

which has the solution Et{z) = .4 + Bz^. According the dictionary of the correspondence, 
A should be zero in order to study the quasinormal states. Close to the horizon (z 1~), 
the differential equation is given by 

E^z) + ^^E'Az) + Y^^^.Eriz) = , (4.3) 

with the solution Et{z) = ^'(1 — z)**^/^ + i3'(l — z)~^^/^. The infalHng boundary condition 
is fulfilled by choosing A! = 0. 

At this point we can perform the next transformation Et{z) = {l — z)~^'^/'^y{z), in order 
to split the infalling singular part from the regular part of the function. In consequence 
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the function y{z) must satisfy the boundary conditions y(0) = and y(l) = 1 and the 
differential equation turns out to be 



y'tiz) + [ai + iuj'yi\y't{z) + [oq + iuj(3i + u;'^ f32\yt{z) = . 



with Qi = Ai, 71 



2(1-2)' 



-ef{z)B{z),p,= '\f;jl-f^ and/?: 



1 



16(1-2) 



(4.4) 
7 + 5(^). 



Results for the quasinormal modes of the transverse vectors in this case are shown in 
figure 4. 

4.2 Longitudinal vectors 

The equation of motion for longitudinal vectors is given by 
■Ai{z){u;^-f{z)k'') + Co{z) 



E'liz) + 



f{z)k^ 



E'l{z) + B{z){u;' - f{z)k')EL{z) = 0, (4.5) 



with Co(-z) = k'^f{z). The asymptotic behavior of this equation is the same for the trans- 
verse e.o.m, then if we do the same transformation that above, we obtain this equation: 
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In the case with k = the differential equations for transverse and longitudinal fluc- 
tuations are the same, in consequence their quasinormal spectra coincide, see figure 4. 

Results for vectors Figure 4 shows the first and second QNM in the complex frequency 
plane. Starting with zero quark mass, i.e. at high temperature (red color), the imaginary 
part monotonously decreases with decreasing temperature. This means the corresponding 
mode becomes more and more stable. In contrast to that the real part of the quasinormal 
frequency first grows until it reaches a maximum and then decreases as well with decreasing 
temperature. This maximum in the real part of the QNM lies above the meson melting 
transition (indicated by a short horizontal dash). The melting transition takes place at a 
critical angle xo = 0.939 from which on the Minkowski embeddings are thermodynamically 
favored, not the black hole embeddings. We have chosen to remain in the so-called under- 
cooled phase keeping the black hole embeddings even beyond the transition. This phase 
is accessible since the meson-melting is a first order transition. So the undercooled phase 
is metastable. However we will see in the following section that at a smaller temperature 
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Figure 4: Location of tiie first (left) and second (right) quasinormal modes in the complex 
frequency plane for the vector fluctuations at vanishing momentum (fc = 0) as a function of the 
embedding xo- Red color indicates small quark mass, or high temperature, while the temperature 
decreases towards blue colors. The horizontal (black) dash indicates the frequency at the first order 
phase transition where the angle is xo = 0.939. The vertical (red) dash indicates the frequecy at 
which the embeddings become locally unstable at xo — 0.962. The modes are followed down to 
embeddings with xo — 0.999875. 
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2.1075 
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2nd QNM 













4.0054 


-3.9976 
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-4.0004 


0.48 


4.0417 


-3.3366 


3.9324 


-3.3386 


0.92 


3.4397 


-1.3093 


3.4397 


-1.3093 



Table 1: Exemplary values for the first and second vector QNM frequencies at fc = for different 
values of xo parametrizing the D7-embedding. See figure 2 for the relation between xo and the 
quark mass Mq. The first pair of values in each row is obtained from the relaxation method, the 
second pair stems from requiring the shooting solution to vanish at the AdS boundary. We find a 
remarkable agreement. 

below the melting transition, i.e. a larger angle Xo^'^^^°^ ~ 0.962, this undercooled phase is 
destabilized by the scalar fluctuation becoming tachyonic. In the figures this is indecated 
by a red vertical dash. 

Exemplary numerical values for the vector QNM frequencies at /c = are given in 
table 1. We find a remarkable agreement between the values obtained with two different 
methods: the relaxation method and the shooting method. The results are in good agree- 
ment for all parameter regions and in all the cases we treat in this work. Therefore we 
exclusively show results produced with the relaxation method from now on. 
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Figure 5: Shoot and relax: Comparison of the shooting method result (green squares) with the 
relaxation method results (black circles) for the location of the first transverse vector quasinormal 
mode at vanishing momentum fc = 0, density d = 0. Along the curves the temperature is varied. 
Our two methods agree very well. 



4.3 Scalar 

In the DBI-action (3.9) we let the 0-angle fluctuate and we split this fluctuation 5Q into a 
product of its singular and regular parts 6&{z) = (1 — z)~*'^/^zy(z). With this change the 
infalling boundary condition at the horizon is translated into y{l) = 1^ and the Dirichlet 
condition at the boundary implies y{0) = 0. The equation of motion for scalar fluctuations 
then reads 

y"{z) + [ai{z) + ici{z)uj] y'{z) + [ao{z) + Mz)uj + b2{z)uj^] y{z), (4.7) 

with 

ai{z) = Ai{z) + ^ , ao{z) = ^ + ^o(^) - B{zfk^f{z) , 



ci{z) = w7-, IT, b2{z) = B{z) 



2(1 -z) ' W 16(1 -.2)2 

{l-Ai{z){l-z))z-2 



4z{l - zf 



Results for the scalar The first and second scalar QNM at vanishing momentum can be 
found in figure 6. The basic behavior is similar to that of the vector modes. Increasing the 
quark mass from zero the real part of the QNM frequency again shows a turning behavior 
moving first to larger values, then to smaller values of Re(w). However, in contrast to the 
vectors, the scalar QNM frequency also shows a turning behavior in the imaginary part 
Im(u;). This means that increasing the quark mass, i.e. decreasing the temperature, the 
corresponding modes first decay faster, then beyond the turning point they decay slower and 

^Notice that the equation for yi^z) is linear and that we can scale y by and arbitrary constant, therefore 
we can always choose the boundary condition at the horizon to be y = 1. 
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Figure 6: Location of the first and second quasinormal modes in the complex frequency plane 
for the scalar fluctuations at vanishing momentum (A; = 0) as a function of the embedding xo- R^d 
color indicates small quark mass, or high temperature, while the temperature decreases towards 
blue colors. The horizontal (black) dash indicates the frequency at the phase transition where 
the angle is xo = 0.939 whereas teh vertical (red) dash indicates the onset of the instability at 



slower as the mode approaches the real axis. Moreover the scalar QNMs do not asymptote 
to the real axis as fast as the vector QNMs do. Instead the scalar QNM frequencies even 
at large masses still have a considerable imaginary part of roughly 1/2. All the values for 
the scalar modes are excellent agreement with the ones obtained previously in [13]. 

The short dash in the figures again shows the location of the known meson melting 
transition where the initial angle is xo = 0.939. Increasing the mass further while staying in 
the black hole phase, we observe a scalar QNM to become tachyonic. This point is marked 
by a vertcal dash in the figures. Figure 7 shows the appearance of an unstable mode 
explicitely. This particular mode is special since it has vanishing real part but it starts 
with an extra-ordinary large imaginary part of the QNM frequency at zero quark mass. 
Increasing the quark mass this purely damped mode moves closer to the real axis until it 
crosses to become unstable at xo = 0.962 corresponding to the maximal mass for black 
hole embeddings of m ~ 1.31. This particular mode has not been observed in previous 
studies because at vanishing quark mass it is located very deep in the complex frequency 
plane near Imcj ~ —8, while for example the first scalar QNM has Imw = —2 at vanishing 
quark mass. In principle there could be an infinite tower of such purely imaginary modes, 
each crossing the real axis at the larger and larger quark mass. However, the accuracy 
of our numerics proved insufficient to establish additional modes beyond this lowest one. 
In any case once this mode has crossed the real axis the D7 brane embedding is locally 
unstable and can not be taken as a (metastable) ground state. This raises the question of 
what is the true ground state in this regime. It might be that there is another type of D7- 
brane embedding that is reached somehow by condensation of the scalar mode. Another 
possibility is that there is simply no locally stable embedding beyond that point. 

4.4 Schrodinger potential analysis 

In this section we present a qualitative analysis of the quasinormal spectrum using the 



Xo = 0.962. 
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Figure 7: Left: The plot shows a purely imaginary quasinormal mode at fc = as a function 
of the embedding. Right: Zooming into the region where the scalar mode crosses the real axis 
becoming tachyonic approximately at xo = 0.96221 

fact that the equations of motion for the fluctuations can be rewritten in the form of the 
Schrodinger equation (see Appendix C for more details) 



where R* is a tortoise-like coordinate. The Schrodinger potential Vs determines the energy 
spectrum E which is related to the quasinormal spectrum by E = uP'. 

At zero baryon density and zero momentum, the potentials for the vector and scalar 
modes are already discussed in [25] and [14], respectively. Nevertheless we include the 
discussion here for completeness. In fig. 8 we present the Schrodinger potential for the 
vector and scalar fluctuations at different quark masses parametrized by xo- In these plots 
we observe an infinite wall in the potential at R* = which corresponds to the AdS 
boundary. In addition to this wall the potential for the vector modes develops a step- 
shape as we increase the quark mass. In [25] it is shown that the imaginary part of the 
quasinormal frequency decreases as the step gets longer which is consistent with our result 
found in fig. 4. 

For the scalar modes a negative well arises in the Schrodinger potential. This well 
becomes deeper and wider as we increase the quark mass and therefore support a 'bound' 
state with E <Q which corresponds to a tachyonic quasinormal frequency Im w > 0. This 
well and the 'bound' state are studied in [14]. The Schrodinger analysis clearly shows the 
existence of a tachyonic mode which we already found in fig. 7. 

4.5 Discussion: Tachyon and De-singularization 

In this section we discuss why the so-called undercooled phase shows unphysical meson 
spectra which do not approach the known ones in the supersymmetric limit. Further we 
discuss that finite density cures this behavior by de-singularizing the geometry, i. e. by 
smoothing out the limiting embedding. 

As mentioned above the scalar fluctuation becomes tachyonic once the quark mass 
parameter has reached its maximum as a function of xo- It is not to be expected that 
the region beyond that point contains physically relevant or meaningful signatures. This 
region contains the limiting embedding which only touches the horizon and geometrically 



+ Vs^l^ = Eip 
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separates Minkowski from black hole embeddings. Here meson spectra had been studied 
earlier [14,25]. These meson spectra display a singular behavior in the sense that all the 
quasinormal modes (first, second, . . . ) approach one single 'attractor' energy (or frequency) 
near the limiting embedding. The geometric reason for this is the scaling symmetry for 
the embeddings in the near-critical region [12]. That scaling symmetry implies that near 
the critical quark mass (or temperature) there exists no preferred scale on the brane. In 
this sense there is no scale which could determine the distance between resonances in the 
brane fluctuations, i.e. between the distinct meson mass resonances, or quasinormal modes 
equivalently. 

At finite fixed baryon density however, this particular scaling symmetry is broken [16]. 
Therefore the chemical potential introduced together with that density does set the scale 
for the separation between the distinct quasinormal modes. This is somewhat analogous 
to the behavior of the thermal quasinormal mode spectrum when the temperature of the 
dual black hole background is lowered towards zero. In that case the quasinormal modes 
are known to approach one single value, i.e. the spectrum becomes singular in our sense. 
There the scale which determines the distance between distinct excitations is clearly the 
temperature since it is the only scale in the theory. In the zero temperature limit this scale 
vanishes and the scaling symmetry is restored. 

As observed earlier [30] the meson spectrum at finite density approaches the supersym- 
metric one at low temperatures. In this sense the theory is de-singularized by finite baryon 
density. We will see this explicitly in the quasinormal mode spectra at finite density in 
section 6. 



5. Finite momentum but vanishing density 

We now turn to the case of non-zero spatial momentum. 
5.1 Transverse vectors 

We obtain results at finite momentum from the numerical solution of equation (4.1) with 
distinct no n- zero values of k. 
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Figure 9: First quasinormal mode for transverse vector fluctuations at different spatial momenta 
k. Horizontal short dashes across the curves show the location of the meson melting transition. At 
A: « 4.16 a spiral structure appears and the curves also asymptote to a distinct small temperature 
value from that k on. Along the curves the quark mass parameter xo = cos Oq is varied. See figure 2 
for its relation to the quark mass and temperature. 



Figure 9 shows the first of the transverse vector QNMs at different values of the spatial 
momentum k. The behavior in the region A; = 0, ... ,4 is very similar to the A: = case. 
Although its trajectory in the complex frequency plane becomes more wavy at larger k, the 
first scalar QNM still starts at quite large real and imaginary parts in order to approach 
the real axis and smaller real parts when temperature is decreased. Distinct curves for 
different values of k within numerical accuracy approach a single limiting value ujq at small 
temperatures. It is interesting to note that the turning point in Re(a;) mentioned in the 
previous section for A; = disappears when k reaches values between A: = 1 and k = 2. 

A quite substantial qualitative change appears at A; = 4.16, where the trajectory of 
the mode in the frequency plane develops one loop. Also the trajectories at higher k have 
this looping behavior. At the same time these curves with one loop do asymptote to a 
single small temperature frequency value loi as well. But this limiting value is distinct 
from the limiting value which is approached by the low k curves without the loop, i.e. 
uJi 7^ Wo- This fact suggests that the loop-behavior and the distinct limiting value are 
somehow related. For the longitudinal vector fluctuations (see figure 11) we will explicitly 
see that this relation generalizes to all fluctuations and to higher loops in this way: All 
the first QNM trajectories for scalar and vector fluctuations with k < k^ have n loops and 
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they asymptote to a small temperature limit frequency a;„ (within numerical accuracy) 
with Un+i > ojn- Note that the loop behavior appears in a physical, thermodynamically 
preferred phase, i.e. before the meson melting transition and also way before the tachyon 
appears. 

The loops are absent in the second quasinormal mode as seen from figure 35. Never- 
theless the second quasinormal mode also asymptotes to distinct low temperature limits 
u;„, above distinct certain momenta kn- 




Figure 10: Dispersion relation for the first transverse quasinormal mode at distinct quark masses, 
or equivalently temperatures, parametrized by the embedding parameter at the horizon, xvi- See 
figure 2 for the relation between xoi the temperature and quark mass M^. 



Figure f captures the dispersion relation of the first transverse vector QNM at differ- 
ent values of xo- 

5.2 Longitudinal vectors 

We obtain results at finite momentum from the numerical solution of equation (4.5) with 
distinct no n- zero values of k. 

The behavior of the longitudinal vector QNMs is qualitatively similar to that of the 
transverse ones discussed in section 5.1. The mentioned loops in the frequency plane 
trajectory do appear at smaller values k ~ 3.4 in the longitudinal channel than they do in 
the transverse one. However, while in the transverse vector case the loops appeared before 
the meson melting transition, in the longitudinal case the transition takes place before the 
first loop is terminated as can be seen from the figure 11. Just as for the transverse vectors, 
also the second QNM of the longitudinal vectors does not have any loops in its complex 
frequency plane trajectory, as figure 37 shows. 

The dispersion relations for the first and second QNM of the longitudinal vector fluc- 
tuation are depicted in figure 12 and figure 38 in appendix E, respectively. 

5.3 Scalar 

We obtain results at finite momentum from the numerical solution of equation (4.7) with 
distinct non-zero values of k. 

The general behavior of the scalar QNMs is qualitatively similar to that of the longi- 
tudinal vector QNMs. Figure 13 shows the first of the scalar QNMs at momenta between 
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Figure 11: The first quasinormal mode for longitudinal vector fluctuations at distinct values of 
k. Along the curves the quark mass parameter xo — cos 60 is varied. See figure 2 for its relation 
to the quark mass and temperature. 




Figure 12: Dispersion relation for the first quasinormal longitudinal vector mode fluctuation at 
distinct values of the quark mass parameter xo (see figure 2 for its relation to the quark mass at 
vanishing density). 



k = 1 and k = 5. Also in this case the overall behavior is that the real and imaginary 
parts decrease as temperature is decreased along the curves. At small k, e.g. k = 1, there 
is a turning point present in the real as well as in the imaginary part. These turning 
points again disappear between k = 1 and A; = 3. Just like for the first of the longitudinal 
vector QNMs multiple loops form successively for larger values of k. The meson melting 
transition appears before the first of the loops has terminated. Again the number of loops 
seems to be directly related to the low temperature value a;„ to which the curves for all 
k > kn asymptote. Dispersion relations for the first scalar QNM are shown in figure 14. 
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The corresponding figures for tlie second scalar QNM are figure 39 and 40 in appendix E. 
Just like for the vectors there are also no loops in the second QNM for the scalars. 




Figure 13: First scalar quasinormal mode at distinct momenta k. Along the curves the quark mass 
parameter xo — cos Qq is varied. See figure 2 for its relation to the quark mass and temperature. 




Figure 14: Dispersion relation for the first scalar quasinormal mode at distinct values for the 
quark mass parameter xo (see figure 2 for its relation to the quark mass at vanishing density) . The 
real part of the QNM frequency is shown on the left, the imaginary part on the right. 



5.4 Schrodinger potential analysis 

Just as in the previous section at zero density and momentum, we here compute the 
effective potential for the scalar and vector fluctuation equations (4.7) and (4.5), (4.1) (see 
Appendix C for more details). 

We begin by examining the Schrodinger potential for the scalar fluctuations in figure 15. 
The lowest (red) curve shows the potential at vanishing momentum and density at xo = 
0.9999. That is near the limiting embedding, far beyond the thermodynamic transition to 
Minkowski embeddings and far beyond the appearance of the tachyon in the spectrum. The 
results at smaller xo are qualitatively the same, but we plot large xo in order to investigate 
the tachyon and the reason for the different "attractor" frequencies in that large mass 
regime. The scalar potential clearly exhibits a wide negative dip in which the tachyonic 
scalar mode resides, compare figure 7. As the momentum is increased the potential is 
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Figure 15: The scalar Schrodinger potential Vs (zooming in on the minimum of the po- 
tential) versus the radial coordinate R* defined in section 4.4 at increasing momenta k = 
0, 2, 4, 6, 8, 10, 12, 14, 16 from bottom to top curve with the quark mass parameter xo = 0.9999 
(see figure 2 for its relation to the quark mass at vanishing density), d = Q. The dip supporting the 
tachyon narrows. 

lifted and the negative dip is narrowed. In this way the lowest possible excitation is pushed 
towards more positive energy values becoming non-tachyonic at large momenta. However 
the theory is clearly unstable against condensation of the scalar fluctuations already at 
Xo ^ 0.962 and A; = 0. In figure 16 we zoom out to larger values of the potential. For 
increasing momentum, a step forms near the boundary. When comparing to figure 8 (a) 
we see that the scalar potential at finite momentum is similar to the vector case at zero 
momentum. In figure 16 the step becomes higher and longer for increasing momentum, 
while its plateau becomes shorter, i.e. most of the plateau is located near the boundary. 
Therefore conceptually the analysis of [25] as discussed in section 4.4 applies as in the vector 
case: When the step becomes longer, the imaginary part of the quasinormal frequencies 
decreases. This is consistent with our observations in figure 13 (compare for example the 
initial points where Xo = for increasing momenta from curve to curve). The increasing 
real part of the quasinormal frequency observed in figure 13 is due to the fact that the 
potential step rises closer to the horizon at larger momentum. This narrows the potential 
and rises the excitation energies. 

Turning now to the transverse vector fluctuations, we observe a step potential in fig- 
ure 17. The larger the momentum, the earlier the potential rises towards infinity when 
approaching the boundary at R* = 0. So effectively the boundary moves towards the 
horizon and the length of the plateau of the potential step becomes shorter. This is differ- 
ent from the scalar modes discussed above. While the imaginary part of the quasinormal 
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Figure 16: The scalar Schrodinger potential Vg (zooming out to larger values of the potential) 
versus the radial coordinate i?* at increasing momenta fc = 0, 2, 4, 6, 8, 10, 12, 14, 16, 18, 20, 40 
from bottom to top curve with quark mass parameter xo = 0.9999 (see figure 2 for its relation to 
the quark mass at vanishing density), d = 0. A step forms towards the boundary, similar as that 
for the vectors. The grid line at i?* = 1 serves to guide the eye only. 

frequency decreases for increasing momentum just as in the scalar case, the Schrodinger 
potential shows a different behavior: For increasing momentum, the Schrodinger potential 
approaches the shape of a wall. Where as in the scalar case, the formation of the step is 
responsible for lowering the imaginary part of the quasinormal frequency, here we expect 
that the formation of the wall is responsible for a similar decrease of the imaginary part of 
the quasinormal frequency. 

5.5 Discussion: Breakdown of hydrodynamics, 'attractors' and the tachyon 

In this discussion we focus on three distinct physical implications of the quasinormal modes 
described previously in this section. The longitudinal vector modes tell us when the hydro- 
dynamic approximation breaks down, while the novel purely imaginary scalar mode renders 
the whole theory unstable as it becomes tachyonic. Both scalar and vector modes asymp- 
tote to 'attractor' frequencies. This behavior is probably related to a spiraling behavior of 
the quasinormal mode's trajectories with changing quark mass or temperature. 



5.5.1 Hydrodynamics to Collisionless Crossover 

Contrary to the transverse vector and scalar channel the longitudinal vector channel also 
has a hydrodynamic quasinormal mode, i.e. a mode whose dispersion relation does not 
show a gap at zero momentum, \\m.k^oj{k) = (see figure 18, left). This mode represents 
diffusion of baryon charge. It is a mode whose frequency is purely imaginary and therefore 
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Figure 17: The transverse vector Schrodinger potential Vs versus the radial coordinate at 
increasing momenta fc = 0, 2, 4, 10, 20, 40 from bottom to top curve with quark mass parameter 
Xo = 0.9999 (see figure 2 for its relation to the quark mass at vanishing density), d = 0. The 
potential step becomes shorter since the potential rises farther and farther away from the boundary 
i?* = 0. 

results in a purely damped time evolution without any oscillation. At small momentum 
the dispersion relation is well approximated by the diffusion kernel uj = —iDk"^, as seen 
in figure 18 on the left. Fitting our numerical data to this we can extract the diffusion 
constant D. Since it has been calculated before in the literature in [14] we do not restate 
this result. We checked however that our values are consistent with the results there. 

On general grounds one expects that a many body system shows a crossover from hy- 
drodynamic behavior at long wavelengths to a coherent or collisionless behaviour at small 
wavelengths. In the holographic context this has first been discussed in an AdS^ example 
in [31] by calculating spectral functions. A more direct way to see this crossover can be 
obtained by studying the quasinormal mode spectrum. At small momentum the hydrody- 
namic mode should be the dominant one, i.e. the one with the largest imaginary part. At 
small wavelength we expect the dominant modes to have frequencies whose imaginary part 
is much smaller than their real parts. This means that at large wavelength the response 
would be simply an exponential decay whereas at short wavelength the response would be 
a slowly decaying oscillation. In terms of the quasinormal mode spectrum this implies that 
the purely imaginary diffusion mode as a function of momentum has to cross the imagi- 
nary part of the dispersion relation for the lower non hydrodynamic modes. Indeed this 
is what happens for the R-charge and momentum diffusion in the strongly coupled M = 4 
theory as discussed in [27,28]. We therefore define the crossover from the hydrodynamic 
regime to the collisionless regime through the momentum at which the imaginary part of 
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the lowest non-hydro dynamic mode crosses the purely imaginary diffusion mode. From 
that wavelength on it is the lowest gaped quasinormal mode which dominates the late time 
response^. It should be mentioned that there is at least one other way of how this crossover 
can be established in terms of quasinormal modes. It could also be that the purely imagi- 
nary diffusion mode pairs up with another purely imaginary but non- hydro dynamic mode 
which allows them to develop real parts as well and to move off the imaginary axis. This 
seems to be the preferred mechanism for ^4^54 black holes [33] and it also appears on probe 
D-branes representing defects in the four-dimensional strongly coupled CFT [34]. As will 
be shown in a companion publication this is also what happens if finite baryon density is 
introduced [35]. This crossover has recently also been investigated for AdS black holes in 
various dimension in the time domain in [36]. 

We have numerically determined the crossover point defined above for different em- 
beddings. As seen from the right side of figure 18, the crossover moves to higher momenta 
as the embedding angle xo is increased. This means that the brane responds to baryon 
charge fluctuations in a purely absorptive way for smaller and smaller wavelengths as the 
quark mass is increased, or equivalently as the temperature is decreased. Of course the 
actual rate of absorption given by the absolute value of the imaginary part of the frequen- 
cies decreases with decreasing temperature. Nevertheless it is somewhat surprising that 
the crossover towards the collisionless regime takes place at smaller wavelengths for lower 
temperatures. 

5.5.2 "Attractor" frequencies 

Here we briefly discuss the appearance of spiraling QNM trajectories and their relation to 
the "attractor" frequencies found in section 5. 

First we should note that in most cases the spirals in the QNM trajectories occur before 
the scalar in the spectrum becomes tachyonic. Thus the spiral is a physical signature on the 
stable or metastable branch of the theory. However the "attractor" frequencies to which 
the trajectories asymptote at large quark mass parameters xo ~^ 1 are located deep in the 
unstable phase of the theory. Therefore the "attractor" frequency is no signature of the 
physical stable sector of the theory. Nevertheless it would be interesting to understand the 
apparent direct relation of these "attractor" frequencies and the number of spirals in the 
QNM trajectories since the spirals are physical, as noted before. 

®There is a small puzzle related to that. If the prolongation of the hydrodynamic mode to large momenta 
is constantly increasing, the front velocity computed from it seems to violate causality. As is well-known 
already the diffusion kernel violates causality because of the behaviour and the extension to larger 
momenta shows even higher exponents in the dependence on k. So how manages the theory to preserve 
causality? The resolution has been presented in [27] : the residues of the diffusive quasinormal mode vanish 
for large momenta and therefore this mode ceases to exist in the dangerous limit of large k. On the other 
hand one can study the hydrodynamics by fixing a real frequency and then looking for complex roots 
in the momentum k as it has been done in [32]. The quasinormal frequencies or the complex momenta 
respectively are roots of infinite order polynomials (or at least of extremely high order polynomials in a 
truncated approximation of the holographic Green function). Therefore it is not possible to simple infer 
the behaviour of the complex momentum modes from the quasinormal modes. Indeed as shown in [32] the 
front velocity of hydrodynamic mode as calculated from the complex momentum roots behaves perfectly 
causal and approaches 1 to very good numerical accuracy. 
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Figure 18: Left: The dispersion relation for the diffusion mode at distinct values for the 
mass/temperature parameter xo (see figure 2 for its relation to the quark mass at vanishing den- 
sity). Right: Intersection of the diffusive mode with the imaginary part of tlie first longitudinal 
quasinormal mode. The intersection point moves to larger values of k, but to smaller values of 
Im^Lu) as the mass/temperature parameter xo is increased. 

Note that a spiraling behavior for changing the temperature has been observed in the 
quark condensate in this system for near-limiting brane embeddings in [12]. There the 
spirals are due to oscillations of the parameters of the embedding, i. e. the quark mass 
and the quark condensate. In particular the asymptotic value m oscillates. This behavior 
has only been observed in the near-limiting embeddings. In contrast to that our loops in 
the QNM-trajectories appear way above the critical embedding already. Nevertheless, as 
stated before there is an apparent connection between the number of loops in our QNM- 
trajectories and the near-limiting 'attractor' frequency. In this way we could argue that 
at finite momentum we see the near-limiting embedding oscillations reflected already in 
the non-critical region in spiraling QNM-trajectories. In other words both the QNM-loops 
and the spiraling quark condensate might have the same origin, namely the oscillating 
embedding parameters which are related to the aforementioned scaling symmetry of the 
near-limiting embedding [12]. 

Unfortunately the Schrodinger potentials at finite momentum but vanishing density 
close to Xo = 1 ill figure 15 and 17 do not show any distinct feature hinting neither on 
discrete special frequencies cOn nor on the jumps between them at critical momenta. 

5.5.3 Tachyon: A new hydrodynamic mode 

We briefly discuss here the behavior of the scalar mode becoming tachyonic as explained 
in the previous section. This mode turns into a hydrodynamic mode in a special case. 

As expected this scalar mode becomes tachyonic at higher and higher values for the 
quark mass parameter xo as the momentum of the excitation is increased. In the parame- 
ter space spanned by charge density, temperature and quark mass there is one interesting 
special point: that is the location Xo^**(^) where the scalar mode becomes tachyonic. Just 
at this special quark mass/temperature value, this scalar mode develops a hydrodynamic 
dispersion relation, i.e. lim/j^o — > 0. In other words the scalar mode which shows the 
instability of the system turns into a hydrodynamic mode just at the critical point. This 
could signal that there is a transition to a new phase. For instance this transition might 
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be similar to the glass transition in supercooled liquids discussed e. g. in [37] . Lastly there 
might not exist a new stable ground state since we might be scanning a regime where no 
stable brane embedding exists besides the thermodynamically preferred Minkowski embed- 
ding. 



6. Finite density but vanishing momentum 

In this section we turn back to zero momentum, but switch on a finite baryon density and 
chemical potential. The non-normalizable mode of the zero component of the gauge field 
living on the D7-brane gives the chemical potential in the field theory. Working in the p 
coordinates we define 

lim Aq{p) = IX . (6.1) 

(0— »oo 

In order to study vector mesons in this background we consider fluctuations of the gauge 
field about this background field. The equation of motion for the zero component of the 
gauge field thus reads [16] 

dpAQ = 2d I = . (6.2) 

/(l-X^) [p^P{l-x'r + 8d^^ 

Since the asymptotic behavior for the embedding function x is known we can also give the 
asymptotics for the gauge field at infinity 

where p, is, as already mentioned, the chemical potential and d parametrizes the baryon 
number density by 

d=^^. (6.4) 

We divided this chapter into four sections in the same manner as the former ones, 
except that wc do not have to handle longitudinal vectors here since A; = 0. But we will 
investigate the transverse vectors and the scalar modes. 

Then wc will describe our results, especially the tachyon behavior and the turning 
points of the meson mass. Furthermore we present a Schrodinger potential analysis to 
understand the above results in a different way. We investigate the turning behavior of the 
quasinormal frequencies further by an analytical calculation close to the horizon at large 
frequencies. 

6.1 Transverse vectors 

To obtain the equation of motion for the transverse vector modes, we have to vary the DBI 
action with respect to the gauge field fluctuations about the background and retain these 
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terms up to second order [30]. The equation of motion for the fluctuations of the gauge 
field reads 



= E +a,,o,[^jVT^^^^^^—^^^jE +2. j,^^,^^-^E. (6.5) 

The connection between the gauge invariant field E and the gauge field is simply E = ujphAi, 
where we are free to choose any Minkowski spatial direction i = 1,2,3. Furthermore 
the parameter a; is a dimensionless frequency defined as w = Uph/i^T with the physical 
frequency ujph. A more detailed derivation of this equation of motion can be found in [30]. 

The solution to this equation of motion can be obtained by numerical methods, de- 
scribed in appendix A. Then we can compute the correlator in the complex plane where 
the quasinormal modes appear as poles. We aim for an explanation of the turning behavior 
of the meson mass found in [30] : 

At large quark mass the meson mass increases proportional to the quark mass Mq as ex- 
pected from the formula found in the supersymmetric case for vector and scalar meson 
masses (for vanishing angular momentum i = on the S^) [8] 

Mmeson = 2{n + I) {n + 2) , n G N . (6.6) 



Whereas for small quark mass there exists a region where the meson mass is decreasing 
when the quark mass is raised. 

Results for transverse vectors Figure 19 (a) displays the paths of the first quasinormal 
mode in the complex frequency plane for different densities. The paths are parametrized by 
the quark mass over temperature ratio. When the quark mass is zero, i. e. our embedding 
is flat, the quasinormal mode is located near the point 2(1 — i) and moves towards the real 
axis as the quark mass increases. For small densities we see that the curves turn around 
and tend towards smaller real values for larger quark mass. This behavior disappears when 
we raise the density up to a critical density dc ~ 0.04. 

6.2 Scalar 

In order to compute the spectral function for scalar modes we first calculate the puUback of 
the metric to the D7-branes and expand in the fluctuations of the embedding variables 6Q 
and 5Q corresponds to the scalar and 6(j) to the pseudoscalar excitations. We consider 
time and p dependence of the fluctuations only, since we stay at zero momentum. The 
induced metric then reads 



-2R^ ^ ^dgSedx^dp + R^daSedbdedx''dx^ . (6.7) 
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Figure 19: (a) First quasinormal modes for several densities d. For smaller densities a turning 
point in the real part occurs, which vanishes as the densities rises, (b) The first quasinormal mode 
for vanishing density. The arrows display the residue of the pole at the specific point, particularly 
the direction its phase and the length its absolute value. 



This coincides with the result found in [14] at zero density. 

To obtain a consistent solution at non-zero density, it is necessary to also include 
fluctuations in the gauge field 5Aq since these couple to the fluctuations 5@ of the induced 
metric of the D7-brane [24]. The coupling occurs since both fields transform as scalars under 
the group of rotations SO{3). They cannot be distinguished by the different transformation 
under the U{1) gauge symmetry anymore since this symmetry is broken at finite density. 
We may also think of the embedding scalar as being effectively charged and this explains its 
coupling to the gauge field fluctuations. Then the action in this case differs from the action 
in [14] by the non-vanishing gauge field and gauge field fluctuation terms. The action can 
be found in appendix D. The equation of motion is given by 



= dp[Adp{5e)] + Bu'^ide) + C{5e) 



(6.8) 



where 



A 



B 



C 



2\3 



(1_ ^2 + ^2^/2)3/2 /1 



2\2 



8d2 



p6/3(l_^2)3+g 



3p^//(1-x')a/1 



8d2 



p6j3(l_^2)3+8d2 



1 - X + P X 



,2^/2 



(6.9) 



with / and / defined in equation (3.2). For a more detailed derivation of the Lagrangian 
we refer to appendix D. 
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Figure 20: The scalar quasinormal modes with pure imaginary quasinormal frequencies for dif- 
ferent densities d. For densities smaller than the critical density dc — 0.00315, the mode becomes 
tachyonic in some region of the parameter space. The several colors corresponds to different den- 
sities d=0 (black), 0.001 (blue), 0.002 (red), 0.0025 green, 0.003 (orange), 0.00315 (cyan), 0.0035 
(purple) . 

Now we can, similar to the procedure with the transversal vector modes, split the solu- 
tion in a regulating and a regular part, and compute the asymptotic solution for the latter. 
This allows us then to apply the shooting technique (see A) to compute the quasinormal 
modes in the same manner as above. 

Results for scalars First we study the density dependence of the mode with purely 
imaginary quasinormal frequencies found at zero density in fig. 7. Our numerical results 
are shown in fig. 20. In this figure we observe that the critical parameter Xo'^'*; which 
the instability occurs, i. e. the quasinormal frequency enters the upper half plane, increases 
with the density. Thus we write Xo^^^C'^)* Also we note that at finite densities the modes 
become stable again at some larger value for xo; which we denote by ('^)- increase 
the density further, we obtain a critical density dc = 0.00315 at which the mode is always 
stable. Therefore the system is unstable in the parameter region Xo"*(^) < Xo < (^) 
and d < dc- The numerical values can be obtained from fig. 20. In section 6.5 as well 
as in the introduction, we further discuss this mode and relate the instability found to an 
instability already known in thermodynamics. 

Let us now consider the behavior of the first quasinormal mode if we vary the density 
d. In fig. 21 we present our numerical results. In this figure we observe that for each 
finite density we find that for large enough xo the quasinormal frequencies behaves in a 
similar way. The real part increases while the imaginary part decreases as we increase 
Xo- For smaller values of xo we observe three distinct movements of the quasinormal 
frequencies. For small densities d < 0.1 the quasinormal frequency first follows the line 
of the quasinormal mode corresponding to zero density as we increase the parameter xo- 
After a critical value of xo is reached the quasinormal frequency leaves this line as the real 
part of the quasinormal frequency increases monotonically. For slightly larger densities 
0.11 < d < 0.2, the frequencies also first moves along the line of the quasinormal mode at 
zero densities but in contrast to the case discussed above the real part of the frequencies 
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Figure 21: Dependence of the first scaZar quasinormal frequency on tiie density: In (a) and (b) the 
quasinormal frequencies for small densities d < 0.2 are shown while in (c) densities up to c? = 5 are 
considered. In (a) and (b) the different colors corresponds to distinct densities d = (black), 0.01 
(green), 0.05 (blue), 0.1 (red), 0.11 (orange), 0.125 (brown), 0.15 (magenta), 0.2 (cyan). In (c) the 
color coding is d = (black), 0.2 (green), 0.25 (blue), 0.5 (red), 1 (orange), 2 (brown), 5 (magenta). 
In all three plots the black dot marks the critical value of of the quark mass/temperature parameter 
Xo where the instability occurs. 

now first decreases. For even higher densities d > 0.2, the frequency at xo = strongly 
depends on the densities. We find from fig. 21 (c) that both the real and imaginary part of 
the quasinormal frequencies at xo = increase with d. As we increase xo we find the usual 
behavior, the real part of the frequency increases while the imaginary part decreases. 

6.3 Schrodinger potential analysis 

In this section we present the Schrodinger potential analysis which we introduced in section 
4 at finite baryon density. We use this analysis to explain the qualitative movement of the 
quasinormal frequencies as we change the baryon density. 

For the vector fluctuations we observe in fig. 19 (a) a turning point in the movement 
of the quasinormal frequencies at low baryon density. The real part of the quasinormal 
frequencies first increases and later decreases as we increase the mass parameter m. This 
turning point disappears as the critical density of dc = 0.04 is reached. After a critical 
quark mass is reached the real part of the quasinormal frequencies always increases while 
the imaginary part decreases. 
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Let us now consider the Schrodinger potential which correspond to these quasinormal 
modes in fig. 22. A similar analysis was also done in [23]. In addition to the infinite 
wall at R* = which corresponds to the AdS boundary an additional peak appears in the 
potential as we increase the quark mass parametrized by xo- For small densities (see fig. 22 
(a)) this peak slowly grows out of the step-shape potential already observed in section 4. 
The step-shape potential is also present at zero density and it is known from the analysis 
we presented there that the corresponding quasinormal frequencies show the turning point 
behavior discussed above. Thus the new feature of the finite density setup is the appearance 
of the peak at high quark masses. As the peak grows, 'bound' states with positive energy 
can be formed, i. e. the real part of the corresponding quasinormal frequency is always 
bigger than its imaginary part. Therefore we find quasiparticle excitations whose masses 
increase as we increase the quark mass. 

If we increase the baryon density, the peak already appears at lower quark mass and 
can thus destroy the step-shape potential (see fig. 22 (b) and (c)). For instance a step is 
still observable in the orange and blue curve in fig. 22 (b) while in fig. 22 (c) this step 
is gone. Since we know that the step-shape potential is the reason for the turning point 
potential, we also find in this analysis a critical baryon density at which the turning point 
disappears. This critical density agrees with the value found in fig. 19 (a). For even larger 
densities (see fig. 22 (d)) the peak increases even faster, i.e. the real (imaginary) part of 
the corresponding quasinormal frequencies increase (decrease) even faster. 

Let us now consider the scalar fluctuations at finite density which show several distinct 
features: a tachyonic mode at small densities and three qualitative different movements of 
the quasinormal frequencies. 

From fig. 20 we know that there exists a parameter region where a scalar mode becomes 
tachyonic. Especially interesting is that at finite density the mode is stabilized as the quark 
mass is increased and that there is a critical density d = 0.00315 at which the mode is 
always stable. In fig. 23 we present the Schrodinger potentials of the scalar fluctuations in 
the relevant density region. As in the case of zero density a negative well appears in the 
potential. It also first grows as we increase the quark mass (see fig. 23 (a)). However as 
we increase the quark mass further, the well decreases while a peak forms in the potential. 
After a critical quark mass is reached, this well cannot longer support a 'bound' state with 
negative energy, i. e. the imaginary part of the quasinormal frequency becomes negative 
again. At the critical density d = 0.00315 (see fig. 23), the potential also shows a negative 
well. However this well just reaches a critical depth in order to support a zero energy 
'bound' state which is due to the zero-point energy (cf. the three dimensional potential 
pot know from quantum mechanics). 

Now we investigate the Schrodinger potentials which are relevant for the movement of 
the first quasinormal frequency studied in fig. 21. These potentials are plotted in fig. 24. 
The first observation is that for small densities d < 0.15 (see fig. 24 (a) and (b)) the 
potentials at low quark mass agree with the potential at zero density. Thus also the 
corresponding quasinormal frequencies must agree which we already found in fig. 21. At 
larger densities (see fig. 24 (c)) even the potential at zero quark mass differs from the 
one at the zero density such that the corresponding quasinormal frequencies at zero quark 
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(c) (d) 

Figure 22: Schrodinger potential of the vector fluctuations at finite baryon density (a) d — 0.003, 
(b) d = 0.01, (c) d ~ 0.03 and (d) d — 0.1. The different colors corresponds to distinct quark 
masses parametrized by xo = (black), 0.5 (green), 0.9 (blue), 0.95 (orange), 0.97 (purple), 0.99 
(red). See figure 3 for the relation between the temperature and the quark mass Mq. 




(a) (b) 

Figure 23: Schrodinger potential of the scalar fluctuations at finite baryon density (a) d = 0.002, 
(b) d = 0.00315. The different colors corresponds to distinct quark masses parametrized by xo = 
(black), 0.5 (green), 0.9 (blue), 0.95 (orange), 0.97 (purple), 0.99 (red). See figure 3 for the relation 
between xoi the temperature and the quark mass Mq. 

mass depend on the density which is consistent with the result found in fig. 21 (c). As for 
the vector fluctuations a peak appears in the potential at finite density as we increase the 
quark mass. These peak can again support 'bound' states which correspond to quasiparticle 
excitations. Thus the imaginary part of the quasinormal frequencies decreases while the 
real part increases as we increase the quark mass. This is the overall movement of the 
quasinormal frequencies which is shown in fig. 21. 
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Figure 24: Schrodinger potential of the scalar fluctuations at finite baryon density (a) d = 0.01, 
(b) d — 0.15, (c) d — 0.5. The different colors corresponds to distinct quark masses parametrized 
by Xo = (black), 0.5 (green), 0.9 (blue), 0.95 (orange), 0.97 (purple), 0.99 (red). See figure 3 for 
the relation between xoi the temperature and the quark mass Mq. 



In general we observe that at finite baryon density a peak appears in the Schrodinger 
potential as we increase the quark mass. If the peak is high enough, it separates the 
horizon of the black hole from the AdS boundary. Only a small leak into the black hole is 
possible due to tunneling. The potential thus approaches the one of a Minkowski embedding 
where the brane does not fall into the horizon of the black hole and therefore can support 
stable normal modes which are calculated in [8] in the supersymmetric limit. In fig. 25 we 
explicitly confirm that the Schrodinger potential obtained from the black hole embeddings 
converges to the potential of the supersymmetric Minkowski embedding^ with the same 
mass parameter m as the mass parameter goes to infinity. This convergence explains that 
also the quasinormal frequencies approach this supersymmetric mass spectrum if the quark 
mass is big enough. This behavior was already found for the vector fluctuations in [23,30]. 
We can also understand this convergence if we look at the phase diagram (see fig. 26). 
For large quark mass over temperature ratios the equal density lines approach the shaded 
region where the Minkowski embeddings are preferred. Since the phase transition is third 
or second order for large quark mass over temperature ratios as it is shown in [21], we also 
expect a smooth transition from the spectrum of the quasinormal modes obtained from the 

''The Schrodinger potential of a supersymmetric Minkowski embedding is given by [25] 

Vs = m [l/4 + 3/8 (tan^ (mi? * /V2) + cot^ (mii * /V2)) 1 . 
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Figure 25: Comparison of the Schrodinger potential obtained from black hole embeddings (solid 
lines) and from supersymmetric Minkowski embeddings (dashed lines) at the same quark mass m. 
The different colors correspond to different quark masses parametrized by xo = 0-99 (black), 0.999 
(green), 0.9999 (blue) at a fixed density d — 0.5. See figure 3 for the relation between xo, the 
temperature and the quark mass Mg. 
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Figure 26: Sketch of the phase diagram: The chemical potential /i divided by the quark mass Alq 
is plotted versus the temperature T divided by M = 2Mq/\/X. Two different regions are displayed: 
The shaded region with vanishing baryon density where Minkowski embeddings are preferred and 
the region above the transition line with finite baryon density where the black hole embeddings 
are preferred. Here we work in the second phase. The curves are lines of equal baryon density 
parametrized by d. 

black hole embeddings to the spectrum of the normal modes obtained from the Minkowski 
embeddings. 

Next we would like to understand the appearance of the peak in terms of the D7- 
brane embedding. In [16] it was shown that the finite baryon density on the brane is 
induced by fundamental strings which are stretched from the horizon of the black hole to 
the D7-brane. At large quark masses these strings form a throat close to the horizon. We 
confirm numerically that the end of this throat and the peak in the Schrodinger potential 
are located at the same value of the radial coordinate p (see fig. 27). Thus we can interpret 
this throat as a potential barrier for the fluctuations which becomes bigger as the throat 
becomes deeper. 

6.4 Analytic solution at high frequencies 

Motivated by the numerical solution to the fluctuation equations of motion shown in fig- 




d^4 y 




'^^vd = 0.00315 
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Figure 27: Comparison of the location of the peak in the Schrodinger potential and the throat in 
the embedding of the D7-brane. The end of the throat is located where x first changes its value. 
The different colors corresponds to different quark masses parametrized by xo = (black), 0.5 
(green), 0.9 (blue), 0.95 (orange), 0.97 (purple), 0.99 (red). See figure 3 for the relation between 
Xoj the temperature and the quark mass Mq. The baryon density \s d = 0.25. 
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Figure 28: Imaginary part of the solution to the regular function F versus the proper radial 
coordinate s [38] . The left plot at vanishing baryon density d — Q shows that the initially sinusoidal 
solution is deformed as the mass parameter xo = 0.01, 0.5, 0.8, 0.9 is increased (see figure 3 for 
the relation between xoi the temperature and the quark mass Mq). Furthermore, its amplitude 
decreases while the wave length increases. The right plot shows that introducing a finite baryon 
density d = 0.2 causes the solutions to change their behavior with increasing xq'- While the first 
three curves for xo = 0.01, 0.5, 0.8 show the same qualitative behavior as those in the left plot, the 
blue curve for xo — 0.9 clearly signals a qualitative change with its increased amplitude. Looking 
at the wave lengths in the lower plot we realize that already the green curve (xo = 0.8) shows a 
decreased wave length as well as the blue curve (xo = 0.9). The awkward oscillation pattern near 
s = stems from a second mode being superimposed on the one we are tracking here. 



ure 28, we suspect that this damped oscillating curve near the horizon can be approximated 
by a damped quasi-harmonic oscillator, i.e. we should be able to find an approximate equa- 
tion of motion which is a generalization of the damped harmonic oscillator equation. By 
quasi-harmonic we mean that the oscillator is damped with the damping depending on 
the location of the mode in radial direction. From the observations in figure 29 we have 
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Figure 29: Left: Solutions to tiie damped harmonic oscillator equation (6.14) for increasing 
damping coefficient 7 = 0, 0.1, 0.2, 0.3, 0.9 (black, red, green, blue, orange). Here the eigenfre- 
quency u) is fixed to 1 in analogy to the plots in figure 28. Right: Hi-frequency {lo = 50) numerical 
solution showing that the damping of the mode is strong near the horizon s = but decreases 
exponentially in the bulk towards the boundary s — * 00. 

already concluded that the amplitude changes rapidly near the horizon and change less in 
the bulk until the boundary is reached. Thus it is reasonable to assume that the damping 
of the mode F already is strong near the horizon and a near-horizon approximation can 
capture this effect to certain extent. In this spirit we take the near-horizon limit ~ 1 and 
at the same time the high-frequency limit ^ 1. 

Applying these limits for the flat embedding xo = in the equation of motion (6.5), 
we obtain the simplified equation of motion 

yH" + {-ioj - y) + + 1) i/ = , (6.10) 

where the variable is y = ilj^^q — 1) and the regular function H(y) comes from the 
Ansatz E = {g — 1)^F with the redefinition F = ^"^^^^^^^^ H . This equation of mo- 
tion has the form of Kummers equation^ which is solved by the confluent hypergeometric 
function of first H = iFi [—iix}{l/ \fl -|- l)/2, — icj, ^ — 1] and second kind U . Boundary con- 
ditions rule out the second kind solution which is non-regular at the horizon and therefore 
contradicts the assumptions put into the Ansatz E = {g — 1)^F. Since we are interested 
in how the solution changes with decreasing m, we need to choose xo non-vanishing. Also 
with this complication we still get Kummers equation with changed parameters and the 
analytic solution for F is given by 



iFi 



+ 1 



, -lUJ, 



9 / 7 I *X2\xoA^ 



(6.11) 



-(.-1)1/4 + ^3^ 




Figure 30: Approximate analytic solution compared to the exact solution at w = 70 , d = , xo = 
0.4 (see figure 3 for the relation between xoi the temperature and the quark mass Mq). 

with the near horizon expansion of the embedding function X = Xo + X2[X0; d]{Q — 1)'^ + . . . 
where we determine recursively 

The approximate solution for F is shown in figure 30. Furthermore we can calculate the 
fraction d^E/E appearing in the spectral function near the horizon using this analytic 
solution. The result is displayed in figure 31. This near horizon limit is not the spectral 
function since we would have to evaluate it at the boundary which lies far beyond the 
validity of the near horizon approximation. Nevertheless, in the spirit of [39] and according 
to our initial assumptions that the effect of damping mainly takes place near the horizon 
we further assume that the limit shown in figure 31 already contains the essential features 
of the spectral function. Indeed the fraction shows distinct resonance peaks which move 
to lower frequencies if we increase the mass parameter m. The right picture shows the 
same situation at a finite baryon density d = 1 and we see that the peaks do not move to 
lower frequencies as much as before. Thus also the vanishing of the turning point at large 
densities as observed before is captured by this approximate solution. 

6.5 Discussion: Turning point and Tachyon 

We now discuss the turning point in the spectrum of the vector modes and its dependence 
on the density. For scalars the purely imaginary scalar quasinormal mode reaches into 
the upper half of the complex plane and thus yields a tachyonic excitation, which we also 
discuss here. 
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6.5.1 Turning point 

In this section we discuss why the vector mesons at finite temperature get smaher masses 
as the quark mass is increased. Further we investigate why there is a turning behavior at 
finite density. In the latter case we also find an analytic solution valid at high frequencies. 

Line of argument Let us first summarize our results. The D3/D7-system at finite 
density has two competing geometrical features. One is the formation of a potential barrier 
near the horizon due to the charge located at the horizon. This causes the resonances 
in the spectral function to become more stable and to move to larger frequency when 
the quark mass is increased. The second feature is the length of the brane ^ supporting 
an excitation from the horizon to the boundary. Increasing this length causes the vector 
meson resonances to move to smaller frequency with increasing quark mass. These two 
geometrical features compete at finite density while at zero density the potential barrier is 
absent. Below we will explain this in greater detail and give a field theoretic interpretation. 

Guiding features of numerical solutions A detailed analysis of the numerical solu- 
tions [38] shows some interesting features. Most prominently there is a turning behavior, i.e. 
the location of the vector resonance peaks first move towards lower real frequencies, then 
turn around and move towards higher ones in order to asymptote to the SUSY-spectrum 
(increasing the quark mass to infinity is equivalent to the zero temperature case). This 
turning point for the lowest of the vector resonance peaks is only present at small densities 
d < 1. For larger densities the peaks run towards higher frequencies with increasing quark 
mass from the start, i.e. there is no turning point. As we see from the vector quasinormal 
modes the imaginary part monotonically decreases with increasing quark mass. At vanish- 
ing density the resonances only move to lower frequencies with increasing quark mass in 
the regime of black hole embeddings. 

Therefore we have two things to explain: 1) the turning of the real part of the frequency 
and the decreasing imaginary part and 2) the motion towards smaller frequencies at zero 
density or small densities and quark masses (d < 1, xo < 0.5). The solutions E of the bulk 

*Since the AdS-boundary is infinitely far away from the horizon, the length of the brane has to be 
renormalized by subtracting for example the length of the brane at trivial embedding xo = 0. 
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Figure 32: Schematic outline of the solution Ax{R*) to the vector fluctuation equation of mo- 
tion (6.5). The horizon is located at i?* = oo where the wave function diverges, the AdS-boundary 
is R*bdy = 0. The solid black curve shows the solution at finite temperature, density and quark 
mass, i.e. on a black hole embedding. The dashed blue curve shows the corresponding solution on 
a Minkowski embedding at the same quark mass. The three regimes correspond to the nature of 
the potential shown in figure 25: I) Minkowski-like, II) potential barrier. III) essentially vanishing 
potential. 

fluctuation equation of motion are composed of a singular part and the regular part F. 
From these solutions the spectral functions are computed. The regular part of the solutions 
has interesting properties, in particular: 

1. The proper distance s on the brane 



increases with increasing quark mass parameter xo- Integrating this expression from 
the horizon, to the boundary gives the proper length of the D7-brane. This also 
changes the wavelength of the fluctuations on the brane, as we will explain below. 

2. As we see from figure 28 the form of the numerical solutions written in the proper 
radial coordinate s shows a damped oscillation. Thus we are motivated to map the 
fluctuation equation of motion to a damped harmonic oscillator equation of motion. 

3. At vanishing density d = a larger quark mass parameter xo induces a stronger 
damping of the solutions F{s) (see figure 28, left). 



p 




(6.13) 
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4. At finite density a larger quark mass parameter xo first induces a stronger, then a 
weaker damping of the solutions F(s) (see figure 28, right). 

5. As seen from the high frequency {uj = 50) solution in figure 29 most of the damping 
occurs close to the horizon while the fluctuation propagates towards the boundary 
with exponential damping. 

6. At (i = for X ~^ 1 black hole embeddings asymptote to m = pxlp-^ptdy ^ 1-3; i-e. this 
geometric construction can not produce the T — > limit in the black hole embeddings. 
Furthermore, long before xo 1 one has to jump to Minkowski embeddings at the 
temperature where the instability appears. In other words: at fixed T you can not 
make the quarks heavier than Mg ~ 1.3T. And there is no smooth transition through 
the singularity at the limiting embedding xo = 1 to the Minkowski embeddings. 

Resonances moving to higher frequencies tow^ards the SUSY spectrum At finite 
density there are two effects present which compete. One effect is the elongation of the 
brane. This is discussed below for the case of vanishing density. A specialty at finite 
density is the appearance of a longer and longer spike in the embedding reaching down to 
the horizon, connecting the Minkowski-like bulk part of the black hole embeddings with 
the horizon. 

The second and as it turns out also the more important effect at finite density is the 
formation of a potential barrier near the horizon at R* = oo, see figure 25. This barrier 
effectively cuts off the horizon from the geometry. Only a small part of the wave function 
"leaks" into the region behind the potential barrier where the black hole is located, that is 
region /// in the schematic solution shown in figure 32. In this region /// the Schrodinger 
potential asymptotically vanishes. Therefore the solution shown in figure 32 first drops 
to low values in order to diverge only very close to the horizon (not shown in the figure). 
Region // is the finite radial distance covered by the potential barrier. Here the wave 
function qualitatively drops exponentially. In region / the main part of the solution is 
located between the boundary and the barrier. Here the potential approaches a pot-like 
form more and more similar to the corresponding potential generated by the Minkowski 
embedding as the quark mass parameter is increased. As this happens, the black hole 
fluctuation solution (solid black curve in the schematic plot in figure 32) approach the 
Minkowski solutions (dashed blue curve in the schematic plot in figure 32). 

As we see from figure 27 the potential barrier moves closer to the boundary as the mass 
parameter xo is increased. Thus the potential gets more narrow and the lowest possible 
excitation is raised to a higher energy, i.e. the real part of the corresponding quasinormal 
mode is increasing. Furthermore the barrier becomes higher such that the corresponding 
excitation becomes more stable, i.e. the imaginary part of the quasinormal frequency 
decreases. In other words, less of the wave function leaks over the barrier into the black 
hole. 

Resonances moving to lower frequencies Let us for simplicity work at vanishing 
density first. All our considerations will also apply to the finite density case. The heuristic 
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explanation for the left-movement of resonances has to do with the proper length of the 
D7-brane. At xo = we have a flat embedding which is the minimal length the brane 
can have when measured in s. Increasing the quark mass parameter xo the D7-brane be- 
comes "longer" in the sense that the proper length computed in the coordinate s increases. 
Therefore we can imagine a solution supported on the brane to be "stretched" together 
with the brane, i.e. its effective wavelength increases. Assuming a constant effective speed 
of light, the effective frequency of this solution has to decrease. In this sense xo here acts 
analogously to a damping coefficient 7 in a damped harmonic oscillator. ^ 

We may quantify this intuition examining the solutions F(s) shown in figure 28. In- 
deed, comparing to the solutions (figure 29) of the damped harmonic oscillator equation 

= dfX{t) + 2uJoldtX{t) + u;o'^X{t) , (6.14) 

with damping coefficient 7, we find that increasing xo in figure 28 resembles the effect 
of increasing the damping coefficient 7 in the analogous harmonic oscillator solutions, 
figure 29. 

This leads us to assume that the fluctuation equation of motion on the D7-brane 
for real values of the frequency w G R can effectively be replaced by the equation of 
motion for a damped harmonic oscillator with an effective eigenfrequency We// and with 
the replacements 7 — > 7(^0) > with 7(xo) a monotonous function and loq — > uj{xo = 0). 
By uj{xo = 0) we mean the effective eigenfrequency corresponding to the solution mainly 
influenced by the lowest of the quasinormal modes. 7(xo) can be seen as the effective 
damping coefficient So the effective solution at finite xo reads 

F{s) oc e-^(^o)'^o^e*''^-//(^o) (6.15) 

with the reduced eigenfrequency 

^effixo) = 1 - lixof ■ (6.16) 

Since < 7(xo) ^ 1 and 7(xo) is monotonous, the frequency of the solution which has 
the eigenfrequency ojq at xo is monotonously decreasing with increasing xo- So the quark 
mass parameter xo effectively acts as a damping coefficient. We might even suspect that 
the embedding x{p) acts as a local damping coefficient depending on the radial location. 
In this way the fact that the damping is strongest at the horizon and vanishes towards the 
boundary (see figure 29, right) is consistent with the fact that the embedding x[p) assumes 
its largest value near the horizon and quickly asymptotes to zero towards the boundary. 

^This heuristic picture neglects the fact that the damping of tlie gravity modes in general is a local 
effect, i.e. it depends on the radial AdS-coordinate (through radially-dependent geometry). Tentatively we 
will assume that we can average the damping effects over the radial coordinate and express them globally 
in an effective damping coefficient independent from the radial coordinate. 

^''This identification of 7(xo) is effective in the sense that the actual quantity appearing in the equation 
of motion (6.5) is x(p) which highly depends on the radial location. So we understand the effective damping 
7(Xo) to be a constant which averages damping effects over the whole AdS-radius. 
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Now let us (as a very crude approximation) choose ujq to be the real part of the first 
vector quasinormal mode on the D7-brane at vanishing xo- The quasinormal mode is al- 
ready damped, i.e. it actually has no real eigenfrequency. But let us nevertheless follow our 
recipe and replace the complicated fluctuation equation by the simple damped harmonic 
oscillator with an effective damping 7(xo) ^-nd effective reduced eigenfrequency w(xo) from 
equation (6.16). The decreasing eigenfrequency io{xo) with increasing xo effectively ex- 
plains the left-motion of the resonances in the corresponding spectral functions computed 
from these solutions ImG^ oc dsF{s)/F{s) . Here we have assumed that the lowest 

resonance peak in the spectral function behaves in the same way (moving to lower fre- 
quencies) as the effective eigenfrequency uJeffixo) of the solution F{s) from which it is 
computed. This intuition we get from the fact that in the exact computation both of these 
frequencies are mainly determined by the behavior of the lowest quasinormal mode. 

These considerations shall serve to give an intuition for the behavior of the resonances. 
To be more precise the resonances are actually influenced by all the quasinormal modes, 
i.e. by their location in the complex frequency plane and by their residues. 

6.5.2 Killing the Tachyon 

The scalar tachyon appearing at zero density is stabilized by introducing baryon charge. 
From a critical density dc = 0.00315 on the scalar does not become tachyonic for any value 
of Xo- Thus the finite charge density d cures the instability and stabilizes the black hole 
phase of the D3/D7 system We have described the scalar quasinormal mode signatures 
in detail in section 6.2. The mechanisms explaining this effect are discussed with the help 
of Schrodinger potentials in section 6.3. The negative potential well supporting the tachyon 
is lifted with increasing charge density. As discussed before the appearance of the tachyon 
is at non-zero density connected to the black hole to black hole first order phase transition 
taking place at finite densities < d < 0.00315 between two distinct black hole phases. 
In particular the tachyon appears on the unstable branch of the free energy diagram of 
the phase transition. This branch connects two metastable branches as shown in figure 1. 
At the critical density dc = 0.00315 both the tachyon and the black hole to black hole 
phase transition disappear. It is not clear what is the physical ground state below dc- As 
was argued in [17] the true ground state might be a mixed phase for which the gravity 
description is not known so far. 

7. Conclusions and Outlook 

Our extensive study of the quasinormal modes of the D3/D7-brane system has revealed 
some interesting relations between previously known and unknown phenomena. As a 
main result we have shown that the system is completely stable above the critical den- 
sity dc = 0.00315. That means that the spectrum of scalar excitations does not contain 
any tachyonic mode. Also in this regime the spectrum of mesonic exictations in the field 
theory -corresponding directly to distinct quasinormal modes- is de-singularized at low 

^^Note that there exist Minkowski embeddings with the same chemical potential and the same quark 
mass, but all states in the Minkowski phase have strictly vanishing density d — [16, 17]. 
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temperature or large quark mass through the exphcit breaking of a scahng symmetry near 
the hmiting embedding. I.e. the different meson excitations behave in accordance with 
the mass formula derived for the supersymmetric case [8]. Furthermore for the regime 
below the critical density, i.e. for d < 0.00315 we have established the connection between 
the black hole to black hole phase transition at finite charge density on one hand and the 
tachyonic scalar on the other. Using the Schrodinger formulation of the problem we have 
explained the movement of scalar and vector quasinormal modes in the complex frequency 
plane in great detail. A universal feature of all Schrodinger potentials at finite charge den- 
sity is that they develop a barrier near the black hole horizon which hides the horizon and 
the black hole from the boundary. In consequence the dissipation decreases with increasing 
density and the quasinormal modes asymptote to a normal mode behavior. This behavior 
has also been observed in [40] in a setup where a black hole in AdS^ develops scalar hair, 
i.e. there is also a charge density distributed near the horizon. 

At vanishing baryon density but finite momentum we found a critical wavelength at which 
the hydrodynamic approximation explicitly breaks down. Below that wavelength the sys- 
tem at late times is no longer governed by hydrodynamic but by (propagating) collisionless 
modes only. A few unresolved issues remain for zero baryon density. In that case we found 
a spiraling behavior of the first quasinormal mode's trajectory for both scalars and vec- 
tors. The number n of loops appears to be directly related to what we coined "attractor" 
frequencies uin to which all QNM-trajectories asymptote if their momentum kn lies in a 
certain momentum regime kn-i < kn < kn+i- A direct relation remains unrevealed. In 
the same case we found that a scalar becomes tachyonic at low enough temperature. This 
instability might give rise to a condensation process or more generally to a phase transi- 
tion as suggested by the hydrodynamic behavior of the scalar mode where it just becomes 
tachyonic. It is interesting to ask what the new phase could be and if it exists at all. 
The analysis presented here can straightforwardly be extended to the D3/D7 setup with 
finite isospin density or with spontaneously broken symmetry [41,42]. In these setups there 
is a variety of modes to be studied, in the latter most prominently the fluctuations of the 
order parameter phase which correspond to Goldstone modes becoming new hydrodynamic 
modes in the phase with spontaneously broken symmetry. 
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A. Shooting method 



We now discuss an improved shooting method for solving the equation of motion for the 
fluctuations at complex frequencies. The problem arises when using the standard method 
for solving the equations of motion by just numerically integrating the differential equation. 

As an example for the failure of this naive method we look at figure 33, which shows the 
spectral function for the transverse vector modes for vanishing quark mass and vanishing 
density. Note here that we use the dimensionless frequency Id = Uph/i^nT) = uj/2 and 
momentum t = kph/{'2nT) = k/2. The line of poles found at Im tr = —1 is definitely 
not correct since in [43] the pole-structure of the very same configuration was analytically 
determined to be 

ft) = n(l - i) for n G No . (A.l) 

It turns out that the numerical errors hide the quasi-normal modes behind this „waU" 
at Im ft) = — 1. This wall could be misinterpreted as a branch cut. 

The problem is, that the equations of motion we investigated are not regular at the 
horizon where the initial conditions are imposed. One therefore splits the solution in a 
regular and regulating part: E[p) = {p — F{p) by computing its Frobenius expansion 
near the horizon. Nevertheless, we have to move the starting point for the numerical 
integration slightly away from the real horizon to, say p = 1.00001 where we choose this 
value so close to 1 that we can expect only small deviations from the exact solution. 
However, for values of Im tt) < — 1 the boundary condition yields 

(p-l)--- =(1x10-5)'°^"'-'^^". (A.2) 

This is problematic when Im ft) is smaller than —1, the initial condition is becoming 
large, leading to round off errors. These change the value of the fluctuation at infinity 
dramatically and are responsible for the invalid results of the spectral function in the 
specific region. 

The cure is to shift the starting point for the numerical integration away from the 
horizon, and thus dealing with not so small values in the basis of (A.2), resulting in better 
numerical initial conditions. Of course, the error from starting the integration further 
away from the horizon has to be compensated. This is simply done, by calculating the 
asymptotics of the gauge field fluctuations in the vector case, or the embedding deviations 
in the scalar case at the horizon to higher order. 

When reimplementing the numerical integration with a starting value of e.g. p = 1.1 
and the new asymptotics, the numerical aberration is prevented successfully and one can 
gain insight farther into the complex w-plane. 

A look at the results shows what can be achieved with this more sophisticated method. 
A comparison between the surface plot of the old method (figure 33) and the improved one 
(figure 34) reveals that the expected pole structure can now be seen clearly. In fact it can 
be checked that the position of the poles agrees with the analytical result (A.l) very well. 
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Figure 33: Breakdown of the standard numerical techniques at Im ir ^ — 1. A series of spikes can 
be seen, being a firmly erroneous solution in view of the analytic result (A.l). This problem will 
be solved by means of our improved method, discussed in the text. 



Figure 34: With the improved method nearly four poles can be resolved for vanishing density 
d = and massless quarks m = 0. For the computation of the initial values the expansion in the 
fluctuation equation were evaluated up to eleventh order. As can be seen, the location of the poles 
fits the analytical solution (A.l) very well. 

B. Relaxation method 

The determination of quasinormal frequencies can be understood as a particular case of a 
two-point boundary value problem. A common numerical method for solving such problems 
is the shooting method, where one solves the differential equations with varying boundary 
conditions at one boundary and searches for a solution that approximates the wanted 
boundary condition at the other boundary within some numerical error. We describe that 
method in appendix A. 

A different approach is provided by the relaxation method which allows to fix the 
correct boundary values on both boundaries. Since this method is less frequently used we 
give a short outline below^^. 

The method is based on replacing the differential equations by a system of finite- 
difference equations (FDE) on a discrete grid. Starting from an ansatz solution obeying 
the correct boundary conditions, one varies the value of the dependent variables at each 
point relaxing to the configuration which provides an approximate solution for the FDE 

^■^A more detailed exposition can be found in [44]. 
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within some given numerical error. In our case, we convert our second order complex 
ODE into a set of four first order equations by separating real and imaginary parts of the 
dependent variables. 

More generally we can consider a set of first order ODE's 

— =5i(a;,yi,..-,y^;A) , (B.l) 

where each dependent variable yi{x) depends on the others and itself, on the independent 
variable x and possibly on additional parameters, like A above. In our case this (complex) 
parameter is the quasinormal frequency. These extra parameters can be embedded into 
the problem by writing trivial differential equations for them 



A 



•3 ' 



— - — - = , smce it IS constant . 
V dx 

We assume that there are n additional parameters and include them from now on into the 
set of dependent variables yi, . . . , ^at, N = N + n. 

The solution to the problem involves N x M values, for the N dependent variables in a 
grid of M points. We also have to fix boundary conditions for the dependent variables. 

The system is discretized as usual 

X ^ ^{xk + Xk^i) , y ^ ^{yk + Vk-i) , (B.3) 

for points in the bulk. One may arrange the whole set of y^'s in a column vector = 
(yi, . . . , yN)k^, where the subscript k refers to evaluation at the point Xk, k = 1, . . . , M. 
With this matrix notation, the system (B.l) can be written as 

= Efc = yfc - yfc-i - (xfc - Xk^i)gk{xk,Xk-i,yk,yk-i) , k = 2,...,M , (B.4) 

where the E^ are the aforementioned FDE's. These are the equations that we need to fulfil. 
Notice there are N equations at M — 1 points, so the remaining N equations are supplied 
by the boundary conditions. We will set ni of them on the left at xi, called Ei, and the 
rest n2 = N — ni at xm, called Ej\/_|_i. 

Now we take a trial solution that nearly solves the FDE's E^. By shifting each solution 
y/c ^ yfc + -^yfe and Taylor expanding in the shift, one obtains the relation 

^ dE ^ 9E 

= Efc(y + Ay) ~ Efc(yfc, y^^i) + V ^ Ay„,fe_i + V — ^ Ay„,fc , (B.5) 

N 2N 

^ -Ej^k = (^Sj^n'^yn,k-1^ + ^ (^Sj^n^yn-N,k^ ■ (B.6) 

n=l n=N+l 

This allows to find the Ayk that improve the solution. First we merge the two differentials 

dEj k „ ^-^j k 

Sj,n = 7^ ■ , Sj n+N = — , U = 1, . . . , N , (B.7) 

dyn,k-i oyn,k 
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in a iV X 2N matrix, for each bulk grid position Xk- For the boundaries the expressions 
follow equally 

-Ej,i = Sj,n Ay„,i = ^yn,i , j = n2 + l,...,N , (B.8) 

n=l n=l y^'^ 



N N 



n=l n=l 



-Ej M+1 = Sj^n ^yn,M = Ay„ Af , j = 1, . . . , 722 , (B.9) 

^ ^ dyn,M 



where n runs in both from 1 to N. The whole {NM x NM) matrix S possesses a block 
diagonal structure. In fact since for reasonable systems N « M S \s a sparse matrix. 
This allows the usage of computer packages in which the solution of linear systems with 
sparse matrices can be done in a very efficient way. The actual solution can now be found 
by an iterative process until a desired accuracy is achieved. As measure for the discrepancy 
of an approximate solution to the actual solution we used 



err 



^ M N 

MN EE 

k=i j=i 



Aybm 



scalevar[j] 



< conv , (B.IO) 



where scalevar[j] is an associated scale for each of the dependent variables (e.g. the value 
at the midpoint or so). The idea is that when that averaged value of the shift to get a 
better solution is smaller than conv, we accept the former values we had as the actual 
solution. In our computations we set conv = 10~^. 

In order to obtain the correct boundary conditions we have used the z coordinate 
system where the horizon lies at z = 1 and the boundary at 2: = 0. We split off the ingoing 
boundary condition on the horizon according to ^{z) = (1 — z)^'^^/'^y{z) . Demanding 
y(l) = 1 and y(0) = gives four real boundary conditions. Since in total we have however 
six dependent variables, counting also the real and imaginary part of the quasinormal 
frequency we need two more boundary conditions. We found it convenient to expand the 
function y{z) in a Taylor series at the horizon and compute also y'(l) which provides the 
additional two real boundary conditions. The typical gridsize we used was consisted of 5000 
points. Finally we note that we have implemented to outlined algorithm in GNU/Octave 
[45]. 



C. Schrodinger potentials 

This appendix shows how to compute the effective potentials, i.e. the Schrodinger poten- 
tials for the scalar and vector fluctuations on the D7-probe-brane. In order to do so the 
linearized fluctuation equations of motion have to be rewritten in terms of a new radial 
coordinate R*. This procedure has been described before (e.g. in [14], [25], ...), never- 
theless we include it here for completeness. For convenience we stick to the notation made 
use of in [14] , and we compute all Schrodinger potentials in the p-coordinates introduced 
in section 3. 

We are not interested in the higher angular excitations on the S^, so we separate the 
fluctuations according to (f){p, S^) = y{p)^{S'^). Let us consider exclusively fluctuations 
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without angular momentum on the S^, i.e. $(5''^) = 1. All the linearized vector and scalar 
fluctuation equations of motion are second order ordinary differential equations and can 
be rewritten in the form 

-^^dp [Hid,y{p)] + [i'H2 + He] y{p) = x^^p), (C.l) 

where Hq, Hi, H2 and Hq are in general functions of p and depend on the particular field 
fluctuation considered. Hg only appears in the scalar fluctuations. For transverse vector 
fluctuations we have 

QPP 

Ho = Hi = ^gPPg--, (C.2) 

gtt 

H2 = ~, He = 0. (C.3) 

gtt 

For longitudinal vector fluctuations we have 

gPP 



Ho = Hi = ^ggPPg^^ (C.4) 

gtt 

H2 = ~, He = 0. (C.5) 

gtt 

For scalar fluctuations the equations of motion do not take such a simple general form 
in terms of metric components since the metric itself contains scalar fluctuations. So we 
choose to write explicitly 

° ~ 8/(1 -x^ + P^X'^) ' ~ (i-x'+p'x")'/' V + /^P«(i-x^)^ ' ^^■''> 

^ ^ 3/V"V8(l-x')'[/V'(l-x')^+8<i2(l-6x2)-144d2/^/yXX'(l-x')'] 
^ 8[/3p6(l-x2)3+8j2]' ' 

where we do not include H2 since the above coefficients are computed in the case of vanish- 
ing momentum 6 = but at finite density d. At finite momentum and at vanishing density 
we have H2 = /2//2 in addition to the d — > limits of Hq, Hi and Hg from (C.6). At 
finite momentum and density the scalar fluctuations couple to the vector fluctuations and 
we will not address this complication in this work. 

In any case we can substitute y{p) = hip in equation (C.l) to obtain 



^ = lt)V- (C.7) 



Introducing the new radial coordinate 

R*= I dp/y^Hoip) , (C.8) 



we can rewrite the first term —HqiP" = — 9r=k^^ + Ho'ip' /2. The special choice of /i = 
Ho^l'^ jHx^l'^ eliminates all the terms containing ■0'. Thus the fluctuation equation of 
motion finally assumes Schrodinger form 

-dR,\ + Vsi^ = Ei, (C.9) 
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with the effective energy E = tv^ and the Schrodinger potential given by 



(C.IO) 



From this general formula the specific scalar, longitudinal and transversal vector Schrodinger 
potentials may be obtained by substituting for Hq, Hi, H2, Hg, h using the corresponding 
values from equations (C.6), (C.4) and (C.2), respectively. 

D. Finite density but vanishing momentum 

To abbreviate the longish Lagrangian we introduced the following notation 



1 2 r ,2/2 

a = l- X, o = a + p X , c 



- , so 1 — c 



p^pa^ + Sd? p^pa^ + 8^2 

In terms of these expressions the Lagrangian, expanded to second order, reads 



(D.l) 



C 



-NfTorVhs'-f 



^ Vby/l — c 



i?4 P a' ,^,2 ,15.? a^{a-hc) 



{dt6er + -p'ff 



&3/2(i_c)3/: 



B? . AT .2 



+ \Y.^tx\dM? + p'p/' 



X 



r^6(l -c)3/2 



{5F^o){dp5e) 



2 

- p'F'\^-^^^{8F,,){5e) + \p'ffa^bVl^cY,Q-Sg^\5F,,f 

- p'ffaVbVT^cF'^^g'^^ (-^^^^i^) 

where Co is simply the Lagrangian without any fluctuations but nonvanishing density 

Co = p^ffaVbVl^c (D.3) 
and the boundary terms Ci and £2 are given by 



Ci = do 



p'ffa^/\' 



m 



and £2 = dp 



3 pV/xx'a / cflx2 



(D.4) 



The symbols denoted by Q with upper indices are the components (diagonal part) of the 
inverse background tensor Q = G + lira'F. The most important ones are 



00 



and Gq 



44 



R%{1 - c) 



(D.5) 
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For the off-diagonal part of the inverse background tensor Q^^ = F^^ = —F'^^ holds, because 
of the antisymmetry of the field strength tensor F, and has the value 



rHjVb{l - c) 

Finally, dFij = di6Aj — djSAi is the field strength for the vector fluctuation. 



E. Results for second quasinormal modes 

In this appendix we collect results for the location of the second quasinormal modes (scalar, 
transverse and longitudinal vectors) at finite momentum but vanishing density. 

E.l Transverse vector fluctuations 

Figure 36 shows the dispersion relation of the second transverse vector QNM at different 
values for the mass parameter ^o- 




Figure 35: Second quasinormal mode of transverse vector fluctuations for distinct values of the 
spatial momentum k. 




Figure 36: Dispersion relation for the second transverse vector quasinormal mode. 
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E.2 Longitudinal vector fluctuations 




Figure 38: Dispersion relation for the second longitudinal vector quasinormal mode fluctuations 
at distinct spatial momenta k. 



E.3 Scalar fluctuations 



-lm(tu) 
4.0 t 




Figure 39: Second scalar quasinormal mode at distinct k. 




Figure 40: Dispersion relation for the second scalar quasinormal mode for distinct mass or 
temperature parameter ^o- 
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